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THE SENSITIVITY OF SYNTHETIC APERTURE RADIOMETERS FOR 
REMOTE SENSING APPLICATION FROM SPACE 


by 


David M. Le Vine 


PREFACE 

Aperture synthesis offers a means of realizing the full potential of microwave remote sensing from space by helping 
to overcome the limitations set by antenna size. The result is a potentially lighter, more adaptable structure for applications 
in space. However, because the physical collecting area is reduced, the signal-to-noise ratio is reduced and may adversely 
affect the radiometric sensitivity. 

Sensitivity is an especially critical issue for measurements to be made from low earth orbit because the motion of 
the platform (about 7 km/s) limits the integration time available for forming an image. 

The purpose of this paper is to develop expressions for the sensitivity of remote sensing systems which use aperture 
synthesis. The objective is to develop basic equations general enough to be used to obtain the sensitivity of the several 
variations of aperture synthesis which have been proposed for sensors in space. The conventional microwave imager (a 
scanning total power radiometer) is treated as a special case and the paper concludes with a comparison of three synthetic 
aperture configurations with the conventional imager. 
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I. INTRODUCTION 


Microwave remote sensing from space offers the potential for measuring many of the parameters important for 
understanding the environment of the earth on a global scale (Butler et ah, 1988; Murphy et aL, 1987), Among the 
parameters which can be measured are sea surface temperature (Wilheit and Chang, 1980), ocean salinity (Thomann, 

1976), soil moisture (Wang et al., 1983), and sea ice concentration (Swift et al., 1985). Microwave remote sensing has the 
advantage that it can be done in the presence of cloud cover, permitting measurement of these parameters in regions inac- 
cessible to visible and infrared sensors. Furthermore, because of the strong sensitivity of microwave radiation to the 
presence of water and the ability of microwave radiation to penetrate beneath the surface, unique information is available 
at microwave frequencies of compliment remote sensing at visible and infrared frequencies. 

However, realizing the full potential of passive microwave remote sensing from space requires putting relatively 
large antennas in space. The apertures required for microwave sensors are large compared to those required of visible and 
infrared sensors because of the much longer wavelength in the microwave portion of the spectrum. For example, antenna 
size is the factor limiting the implementation of an L-band radiometer in space to measure soil moisture, and antenna size 
is an important factor limiting development of a microwave sensor to fill the gaps created by clouds in present-day visible 
and infrared soundings of the atmosphere. 

A possible means of overcoming this size limitation is to use aperture synthesis (e.g., Le Vine and Good, 1983; 
Thompson, Moran and Swenson, 1986). This is a technique in which correlation receivers are used to coherently measure 
the product of the signal from pairs of antennas with many different antenna spacings. For distant sources, it can be 
shown that this correlation function is proportional to the Fourier transform of the intensity of the source at a frequency 
which depends on the spacing between the antennas. By making measurements at many different spacings one determines 
the Fourier spectrum of the source; then a map of the source can be obtained after all measurements are complete by in- 
verting the transform. The resolution obtained in this manner is determined by how well the correlation function has been 
measured, not by the size of the antennas used. In principle, one can obtain very high-resolution maps of the source by 
measuring at many different baselines using relatively small antennas. 

This technique has been successfully employed in radio astronomy to obtain very high-resolution maps of radio 
sources in what is called “earth rotation synthesis" (Swenson and Mathur, 1968; Brouw, 1975; Hewish, 1965; Thompson, 
Moran and Swenson, 1986). The Very Large Array in Socorro, New Mexico is an example (Napier et al., 1983). Several 
versions of this technique have also been proposed for remote sensing from space. These include a proposal (Schanda, 
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1976, 1979) in which the baseline is changed by making measurements at different frequencies, a concept (Mel’nik, 1972) 
involving a moving system with a single fixed baseline in which measurements must be made at several different time de- 
lays, and a modification which employs a scanning linear array which is rotated (C. Wiley, private communications). 

More recent variations include a configuration resembling a Mills cross in which multiple beams are formed with antenna 
elements along each arm of the cross (Milman, 1988); a variation resembling Mel nik s configuration (Mel nik, 1972) us- 
ing contiguous parallel beams and matched filtering (Hughes Aircraft Co., private communications) and finally, a hybrid 
real and synthetic aperture system which uses real antennas to obtain spatial resolution in one dimension and aperture syn- 
thesis to obtain resolution in the other dimension (Le Vine et al . , 1989; Swift et ah, 1986). 

Among the most important issues which must be resolved to determine the viability of aperture synthesis for remote 
sensing from space is the sensitivity (AT) which can be achieved with a particular configuration. Sensitivity is an especial- 
ly critical issue for measurements from low earth orbit because the motion of the platform (about 7 km/s) limits the in- 
tegration time available for imaging a particular scene. The problem is especially pernicious because the higher the spatial 
resolution required the less time is spent by the spacecraft over the scene. Other factors being constant, sensitivity de- 
pends on the actual (i.e. physical) collecting area of the antenna system employed. But the goal of aperture synthesis is to 
reduce the physical collecting area needed for a given spatial resolution. The sensitivity lost in this trade can, in many 
applications, be reclaimed by increasing the integration time-bandwidth product. However, this is not always possible for 
sensors in space because platform motion restricts the integration time, and requirements on the field-of-view and the 
problem of RFI restrict the bandwidth (Thompson and D'Addario, 1982). This trade-off between physical aperture and 
sensitivity was the motivation behind the hybrid real-and-synthetic aperture system mentioned above which was conceived 
as a means for achieving the sensitivity needed for soil moisture measurements from space (Le Vine et al., 1989; Ruf et 
al., 1988). 

The purpose of this paper is to develop expressions for the sensitivity of remote sensing systems which use aperture 
synthesis. The objective is to develop basic equations general enough to be used to obtain the sensitivities of the several 
variations which have been proposed for sensors in space. The analysis is done in two parts: the first is a derivation of 
the sensitivity of the detector used to measure the signal from an individual pair of antennas (a coherent correlation 
receiver) and the second is the derivation of the sensitivity of the pixels in the final image. This is illustrated in Figure 1. 
On the left are shown two antennas viewing the scene. In aperture synthesis the antenna output voltages, V 4 , are multi- 
plied together and averaged. The middle panel of Figure 1 shows the schematic of a detector which makes this measure- 
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ment. The output of this detector, y i} , is a point in the Fourier transform space of the scene; and the scene itself is 
obtained by doing an inverse Fourier transform on the set of measurements {y^}. This is done by the processor (far right) 
in Figure 1. Many variations on this basic theme are possible. The variations may differ only in how the baselines are dis- 
tributed in space or may involve subtle variations in how the signal processing is implemented. However, in each of the 
variations mentioned above, the correlator is a common element and the variations may all be cast in the form of Figure 1 
with different processors. Hence, this paper will begin (Section II) with a discussion of the correlation receiver, deriving 
first the signal-to-noise ratio (SNR) and then the temperature sensitivity (AT) of the output. Then (Section III), an expres- 
sion will be derived for the sensitivity in the Final image, assuming that the processor is a discrete Fourier transform on 
the set of correlator outputs. It will be shown that in the most general case, the sensitivity may not be the same at each 
pixel in the image. Finally (Section IV), these results will be used to compare several of the variations of aperture synthe- 
sis mentioned above as potential candidates for remote sensing from space. The sensitivities of these configurations will 
also be compared with that obtained using a scanning, total power radiometer. The later result is obtained as a special 
case of the analysis done above. 

II. THE CORRELATION RECEIVER 

In aperture synthesis the coherent product (amplitude and phase) of the output voltages from a pair of antennas is 
measured for many pairs at different antenna spacings, and then an image is formed by processing this set. Ultimately, 
one is interested in the signal-to-noise ratio and sensitivity of the image, but as pointed out above, it is convenient to do 
the calculations in two steps, beginning with the receiver and then doing the calculations for the image. In this section, 
expressions are derived for the signal-to-noise ratio and the sensitivity of the correlation receiver. 

A. Signal-to-Noise Ratio (SNR) 

A schematic of the ideal correlation receiver is given as part of Figure 1 (middle). The receiver forms the coherent 
product of X ( and X } which represents the voltages at the output terminals of the k-th pair of antennas. As part of the sig- 
nal processing, these voltages would normally be amplified, filtered and mixed down to a convenient IF for further 
processing. The transfer functions, H t and H ( represent the cumulative effect of this processing. N, and Nj represent the 
equivalent noise in the circuit referred to the input and H LP is a low pass filter which represents the averaging (integra- 
tion) which is done after the multiplication. 

The input X i will be written in the form X, = X t 11 + x } where X t ° is the antenna output from a (perhaps hypotheti- 
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cal) constant background against which changes are measured. The voltage Xj is the change relative to this background 
which it is desired to detect. The signal r ( at the input to the receiver (Figure 1) is the sum of the antenna output and the 
system noise: r = Xj° + N; + x r 

Define the signal-to-noise ratio to be (e.g., Tiuri, 1964) the mean change in output due to the presence of the signal 
divided by the standard deviation of the output when the signal is zero (i.e., divided by the RMS noise). Using this defi- 
nition, the SNR for the receiver output, y ;j (r), is: 


< yTX 0 + N + x) > - < yij (X° + N) > 

SNR \ 0) 

^ < y 2 j j(X° + N) > - < yjj(X° + N) > 2 


Expressions for the moments of y y needed to evaluate Equation 1 can be obtained from Figure 1 using standard 
procedures (e.g., Appendix A). Since the signals r,(t) are the result of thermal noise, it is reasonable to assume that they 
are very broadband compared to the (bandpass) filters Hj(v). In addition, assuming that H^ p (v) is very narrow band com- 
pared to the Hj(v), it can be shown (Appendix A) that: 


< y.jCr) > 


oo 

! 


Sjj(v) H lp (0) \ H,(v) Hj*(v) dv 


Varfy^r)] 


= < yii 2 ( r ) > - < y,,(r) > 2 


= [S,, 2 (v 0 ) + S; 


CO OO 

ii(v 0 ) Sjj(v H ) ] | H 2 (v) H J 2 (v) dv j H lp 2 (v) 


dv 


(2a) 


(2b) 


where S ;j (v) is the Fourier transform of the correlation function, R^( r ) = < r(t) ^*(1 + r ) > , and v 0 is the center fre- 
quency of the bandpass filters. The preceding expressions ignore “fringe washing,” a problem which becomes important 
as the passband of H,(v) becomes large (Thompson and D’Addario, 1982). In employing Equations 2 it is convenient to 
make the following definitions: 
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2 


H,(v) Hj*(v) dv 


[ H,(v) R(v) ] 2 dv 


H 2 L p(0) 


H 2 , p(v) dv 


where B is called the “noise bandwidth” of the receiver (Tiuri, 1964; Thompson et al., 1986) and t_ is the equivalent in- 
tegration time of the low pass filter (Kraus, 1966; Tiuri, 1964). _ 

Now, substituting Equations 2 into Equation 1 and using the definitions in Equations 3, one obtains the following 
expression for the signal-to-noise ratio of the receiver: 


S N R = 


/Sij 2 (Vo) T S„(v o ) Sjj(v 0 ) 


•V 2 B T 


In order to proceed further, the spectral densities, S (j (v), are needed. These are needed in two cases: a) when the in- 
put is the ambient background plus system noise (r f = X,° + NJ, and b) when the input is the ambient background plus 
system noise plus the signal to be measured (r = Xj° + N ( + x It will be assumed that the noise Nj in the two chan- 
nels of the radiometer are independent, and are independent of the ambient noise, X". Then, in the first case (r, = X” 4- 
N|), one obtains the following correlation functions: 

R tj (r>) = < [ X” + NJ [ X/> + N } ] > (5a) 

= < Xj° Xf > 

R n (r) = < [ X 5 ° Hh NJ 2 > (5b) 

= < X,° Xj° > + < N, 2 > 

where < x { x- > has been used as a shorthand notation for the correlation function <x j (t)x j *(t + r)> and the complex 
conjugate (*) has been dropped because the signals are real. If the noise and ambient signal in the two channels are the 
same, the results simplify and one may write: 
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Rii(r') = R/t) = < X G 2 > + < N 0 2 > (6a) 

where X" = Xj° = X 0 and Nj = Nj = N 0 . In the second case = X° + N, + x,), only the cross terms (i * j) are 


needed in Equation 4. One obtains: 

Ry(T) = < [Xj 0 + Ni + x.lfXj 0 + Nj + Xj] > (6b) 

= < X° Xj° > + < Xj Xj > 
where it has been assumed that < Xj Xj° > =0. 

Now, the spectral densities S^v) needed in Equation 4 are obtained by Fourier transforming the R ~(t) in Equations 
6a and 6b. Denoting this Fourier transform by double brackets, < < > >, Equation 4 becomes: 


<<x,x i >> 

S N R = ■ ■ 1 V 2 B r (7) 

< X°Xj 0 > > + [<<X 0 2 >> + <<N 0 2 >>] 2 

The Fourier transforms needed in Equation 7 can be derived following straightforward procedures (e.g.. Appendix 
B; also see Fomalont and Wright, 1974). It turns out that the transforms are proportional to the spatial Fourier transform 
of the scene evaluated at a frequency n = v 0 L/c where L is the distance between antennas. Since Xj° is the output due to 
the ambient background, which by assumption is constant and the same for each antenna, the spectrum, < < X/’Xj 0 > > 
is zero except possibly at 17 =0. Thus, when r\ & 0, the SNR is: 

< < XjXj > > 

S N R = ^2 B r ( 8 ) 

< < X c 2 > > + < < N c 2 > > 

When 17 = 0 the correlation radiometer (Figure 1) is degenerate and it becomes a total power radiometer. In this case, the 
analysis follows as above except that Tj = r and Hj(v) = Hj(v). In this case, one obtains: 


< < x> > > 


SNR 


■ V B r 


< < X , 2 > > + < < N 2 > > 


which differs from Equation 8 only by a factor of 2 (Tiuri, 1964; Kraus, 1966). 


B. Sensitivity 

The sensitivity of the correlation receiver is obtained by setting the left-hand side of Equation 8 to unity (i.e., SNR 
= 1) and solving for the fluctuation < < x,Xj > > required to produce this SNR. However, before doing this it is 
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convenient to change the notation in the preceding expressions to that commonly used in radio astronomy and remote 
sensing applications. 

To proceed, notice that the Fourier transform, < < x i x j > > , can be written in the following form (Appendix B; 
also see Swenson and Mathur, 1968): 


f P „<a> B n^ 

<<x i x l >> = (Z 0 A e /2) 1 / = ex P(— j2?rij»q) dq (10) 

J Vl - q x 2 - q y * ' " - 


where j? = L (v Q /c) and L is the vector distance (i.e. “baseline”) between the antennas, Z 0 = V/i 0 /e 0 is the impedance of 


free space and A e is the effective receiving area of the antenna. (In the general case < < x^ x j > > is a complex num- 


ber; however, its complex conjugate is obtained trivially by interchanging antennas. Hence, one can always form a real 


signal to use in calculating the signal-to-noise ratio by averaging these two values: Re< > = [< <x f Xj> > + 


< <XjXj> >]/2, This value will be used in Equation 8 to obtain the SNR. Additional discussion and an example of this 


point is given in Section 1UB.4.) 


Before substituting into Equation 8, it is convenient to put Equation 10 into a more symmetric form by defining a 


normalized brightness temperature B n (q) as follows: 


B(q) 

B„(q) = = [\-72kT B l B(q) (11) 

B„(q) 

where B 0 (q) is the brightness of a constant thermal source of temperature T B (i.e., assuming that the Rayleigh-Jeans law 


applies to this source). With this definition. Equation 10 becomes: 


C P„(q) B„(q) 

<<X , X J >> = Z 0 [kT B /X 2 )A c \ — — ■ =rexp(-j2Tq«q) dq (12) 

J vi - q , 2 - ~~ ~ 


and defining 6(V) to be the following integral: 


+ 1 


6(yj) = Re 


P„(q) B n (q) 

exp( -j2ir^*q),dq 

Vi - q.x 2 - q y 2 


(13a) 


« 


one obtains: 


Re < < XjXj > > = Z G [kT B /X 2 ] A e 6(v) (13b) 

Notice that when B n (q) = 1 and r\ = 0 (i.e., the special case of a total power radiometer), then d(rfi reduces to the 
“solid angle” of the antenna used in the measurement (Kraus, 1966): 


n 


e 


P n (0, <j>) sin(0) d0 d <f> 

o 


+ 1 



-i 


(14a) 


(14b) 


= X 2 / A e (14c) 

The transformation between the spherical coordinates (0, <j>) and the direction cosines (q x , q y ) is straightforward ( Appen- 
dix D or Swenson and Mathur, 1968) and a derivation of Equation 14c can be found in Kraus (1966). 

It is possible to obtain an explicit expression for < < X 0 2 > >, due to the ambient background, in terms of the no- 
tation introduced above. In particular, since < < X () 2 > > is measured at a = 0 and due to a constant (B n — 1), one has 
< < X 0 2 > > = ZJk T x / X 2 ] A c = Z o k T x . Finally, assuming that the internal system noise is due to an equivalent 
source, but with temperature T N , one can write: < < N 0 2 > > = Z 0 k T N . Now using these expressions and Equations 
13 and 14c in Equation 8, the SNR can be written: 

t b 0(a) # 

SNR V2B r (15) 

t x + t n n e 

Finally, the sensitivity of the correlation radiometer can be obtained by solving for the T B which yields a SNR = I. 


Calling this value AT, one obtains: 



V2 bT 0(5) 


(16) 


where T sys = T x + T N is the total “system” noise. Notice that the sensitivity, AT, depends on the antennas used in this 
measurement (i.e., on their solid angle, fi e ), on the baseline between antennas ( a )» also on the structure of the 
source through 0(2). 
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C. Special Cases 

Equation 16 is a general expression and it will be used below to analyze several potential sensor systems which use 
aperture synthesis. However, before doing so several special cases will be considered to show that Equation 16 includes 
the standard results. 

1. Point Source 

First, consider the case when the scene is a point source. In this case, assuming that the point source is well within 
the main beam of the antenna so that P n (q)/ Vl - q x 2 — q y 2 = 1, one obtains d(ij) =cos(2 i r £*q o ) where q G is a unit 
vector pointing toward the source. Using this result, the sensitivity is: 

T 0 

* sys e 

AT = - point source (17) 

\!l B r cos(0 o ) 

where 0 O = 2 x Notice that the singularity at cos (4> 0 )= 0 is only apparent because cos(0 o )= 1 near the main beam 
of the antenna, where it has been assumed that the point source is located. Equation 17 is identical to that derived by 
other methods (e.g., Tiuri, 1964; Kraus, 1966; chapter 7) except for the effect of the receiving antenna, 0 e , which 
appears explicitly here. This factor doesn’t appear (it is unity) if it is assumed that T sys = T x 4- T N is the noise tempera- 
ture of an equivalent point source (see below) as is done in the references cited above (e.g., Kraus, 1966). 

2, Total Power Radiometer 

As a second check, consider the case of a total power radiometer. The derivation follows as above but begins with 
Equation 9 for the SNR and uses V = 0 (zero baseline) in Equations 10, 12, 13 and 15. The result is: 

T U 

sys c 

AT Total Power (18a) 

VB~r 0 ( 0 ) 

In the special case of a uniform source which fills the beam of the antenna (B = 1), one has 0(0) = fi e and the expres- 
sion for the sensitivity becomes: 

T 

4 sys 

AT — Distributed Source (18b) 

TbT 

which is the classic result (Kraus, 1966; Ulaby et al. , 1981). In this case the sensitivity is independent of the antenna 
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aperture. However, notice that in the case of a point source (located very near the main beam of the antenna) one has 
0(0) = 1 and the sensitivity is: 

T 

AT ft Point Source (18c) 

VbT e 

which depends on the antenna aperture. Clearly, the larger the antenna aperture (ft e = X 2 / A e ) the more sensitive the 
measurement will be to a point source. 

3. Noise Temperature Using Equivalent Point Source 

Finally, note that in each of the examples above, the equation for AT depends on the definition of T B and T N . In 
particular, these could have been chosen to be the temperature of an equivalent point source rather than a distributed 
source (as was done above). Assuming they are due to a point source near the main beam of the antenna (i.e., at — 

0), then < < X 0 2 > > = Z 0 k T x / ft e and < < N 0 2 > > = Z 0 k T N / ft e . The derivation of AT in this case follows as 
above, and the results can be obtained from Equations 15 through 18 by multiplying each equation by l/ft e . For example, 
Equation 17 becomes AT = [T sys /V B r]/ CO s(0 o ) and Equation 18b becomes AT = [T sys />/B r ]/ft e but now Equation 18c 
is AT = T sys /VBT. 

m. SENSITIVITY OF THE IMAGE 

Equation 16 is the sensitivity of the output of a coherent correlation receiver (Figure 1) with an arbitrary, but fixed 
baseline ij. In aperture synthesis, measurements are made at many baselines and then an image of the scene is constructed 
by taking a Fourier transform of these measurements (a discrete Fourier transform on the discrete set of measurements). 

In this section, an expression will be derived for the sensitivity in the final image when the processor in Figure 1 is a dis- 
crete Fourier transform (DFT). This provides a basic formula from which several variations of aperture synthesis can be 
discussed (Section IV to follow). 

For simplicity, assume that the measurements (receiver outputs, y tj ) are made at antenna baselines which can be 
mapped onto the coordinates of a Cartesian grid with uniform spacing, d x and d y along each axis, respectively. That is, 
assume that the measurements can be mapped onto the points ( m d x , n d y ) where m=0, 1,2 — M — I and n = 0, 

1, 2 .... N - 1. Also, let I mn be the intensity in the mn-th pixel of the image. Then, with the processor being considered 
here, I mn is the two-dimensional discrete Fourier transform of the receiver outputs: I mn = DFT [y mn ] - The objective of 
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this section is to derive an expression for the sensitivity of the I mn to changes in the scene. 




A. Signal-to-Noise Ratio (SNR) 

The initial step is to calculate the SNR and then the sensitivity is obtained by setting the SNR to unity. Adopting the 
same definition of signal-to-noise ratio for the image as was used for the receiver, one has: 


S N R =; 


< I m „(X° + N + x) > - < IJX« + N) > 


(19) 


v< I- mn (X<' + N) > - < r mn (X" + N) : > 

In this notation, I mn (X° + N) = DFT [y mn (X° + N)] is the intensity of the mn-th pixel in the image when the input is 
the ambient background (X°) plus system noise ( N ). < I mn (X° 4- N) > is a constant (d.c.) background noise in the im- 
age due to system noise and the ambient scene. The “signal’ ’ is the change in the scene relative to this ambient back- 
ground and is due to x. 

Since the DFT is linear, the numerator in Equation 19 can be found immediately from Equations 2a, 5a and 6b. 

One obtains: 

(20) 


< I mn (X° + N + x) > - < I mn (X° + N) > = DFT{ < < x ni x n > > } H LP (0) 


!"■ 


,(v) H n *(v) dv 


In order to evaluate the denominator more information is needed about the output noise, y mn (X° + N). It is straightfor- 
ward to show that <y mn (X° -F N)> =0 for all non-zero baselines (e.g.. Equations 2a, 5a, and the discussion following 
Equation 7). In addition, it will be assumed that the y mn (X° + N) are all independent and identically distributed. Then 
substituting into the definition of the discrete Fourier transform and using Equation 2b, one obtains: 

< I 2 mn (X° + N) > - < I nm (X° + N) : > = (i/MN) [ <y 2 mn (X° + N)> - <y mn (X° + N)>^ ) (21a) 


= (I'MN) [ S : nm (v il ) - S nn (v„) S min (v tl ) ] 


A 

H m (v) H n *(v) dv 


OO 



H- lp (v) dv 


— oo 


- OO 


(21b) 


where M and N are the number of elements in each dimension of the discrete transform. Notice that the subscript "mn" 
has been kept in the equation above, although by assumption the result is independent of m and n (except possibly at m 
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= 0 and n = 0 where the receiver is degenerate and becomes a total power receiver). Now, using Equations 5a and 5b 
and 6a in Equation 21b and using the definitions given in Equations 3, one obtains: 

DFT {<< x m x n >>} 

S N R = V2MNBT (22) 

< < X 0 2 > > + < < N 2 > > 

Before calculating the sensitivity, it is convenient to rewrite Equation 22 in more conventional notation as was done 
above (Section IIB). Following the same procedure (i.e., defining a normalized brightness B n (q) as in Equation 11 and 
then using Equation 12), one obtains: 


DFT{< <x m x„> >} = Z p [k T b ] [A e /X 2 ] 



G(q mn - q) dq 


= ZJk T b ] [A e /X 2 ]0(m,n) 
where the following definitions have been made: 

G(q mn - q) = DFT {w exp(-j2x^*q) } 


(23) 

(24a) 


f P„(q) B n (q) 

0(m,n) = = G(q mn - q) dq 

J Si - q * 2 - " 


(24b) 


The factor, w, introduced in Equation 24a is an arbitrary “weight” which has been included in the DFT for complete- 
ness. It can be assigned arbitrarily to shape G(q) or can be used to account for redundant measurements. The vector (j m n 
identifies each pixel in the image. Its components q m n = (q m , q n ) in the special case of uniformly spaced baselines are q m 

= m X / d x M and q n = n X / d y N. (See Appendix F). 

Now, assuming that < < X 0 2 > > = Z 0 k T x and < < N 2 > > = Z 0 k T N , as done in Section IIB, and using 
Equation 23b, one obtains the following result for the signal -to-noise ratio: 

Tb 0(m,n) 

S N R = v'2 M N B r (25) 

T x + T n fi e 

where Ii e is the “solid angle” of the antennas used in the measurements as defined in Equations 14. The sensitivity of the 

image is obtained by setting the SNR = 1 and solving for the signal T„ required to achieve SNR = 1. Calling this signal 
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AT, one obtains: 


T sys i n e 

AT = — = — (26) 

\2 B r n/mn 0(m,n) 

Notice that since 0 (m,n) depends on m and n, the SNR and sensitivity, AT, may be different for each pixel. However, 
while this is true in general, it is not a problem for the special cases most commonly encountered in practice. Examples 
will be given in the special cases to be discussed below. 

B. Special Cases 

In the case of uniformly spaced measurements, G(q) has the form of the array factor encountered in the analysis of 
uniform, linear antenna arrays (Appendix F; also, see Jordan and Balmain, 1968). In applications such as remote sensing 
from space, the number of individual antennas in the synthesis array is likely to be large (e.g., Le Vine et al., 1989), in 
which case G(q mn - q) will be a rapidly varying function in comparison to P n (q) with its support near q mn . Hence, near 
the center of the image where | q mn | < < 1 one obtains: 

+ \ 

0(m,n) = | B n (q) n - q) dq (27) 

I 

This is just the expression for the power received by an antenna with the power pattern, G(q), when viewing the source 
B n (q) (e.g., Kraus, 1966; Collin and Zucker, 1969). Hence, close to the center of the image, 0(m,n) will depend on 
(m,n) only to the extent that the scene, B n (q), changes from pixel to pixel. 

Before proceeding, it is convenient to define a solid angle, 0 syn and an equivalent area, A syn , for the “synthesized 
beam.” In light of the discussion above (Equations 24 and 27) and in analogy with Equations 14, the following definitions 
are made: 

p„(q) 

■ G(q) dq (28a) 

1 - q x 2 - q y 2 

X 2 / A syn (28b) 
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With these definitions, Equation 26 can be written in a particularly appealing form in a few important, special cases as 


discussed below. 

1) Distributed Source 

First, consider the case of a distributed source which is constant over the dimension of one pixel. Then 0 (m,n) = 
fi syn = x 2 i A syn and Equation 26 becomes: 


AT 


sys 


J 


syn 


s/2 B T VM N A, 


(29a) 


This is the form commonly found in radio astronomy applications (Napier et al., 1983; Thompson et al., 1986). 


2) Point Source 


Next, consider the case of a point source located in the main beam at q m n . In this case, one has 0 (m,n) n I 


and Equation 26 becomes: 


AT 



V 2 BT 



(29b) 


3) Zero Redundancy 

Additional insight can be obtained by examining the factor, MN, in Equations 29. Consider n antennas and suppose 
that measurements are made at all possible antenna pairs. There are n(n - I)/2 different pairs. If no pairs have the same 
baseline (and the baselines are uniformly spaced), then this arrangement is referred to as a “zero redundancy' 1 array. In 
practice, more than n antennas would be required to achieve the n(n — l)/2 baselines because zero redundancy arrays 
have been found in only a few cases (e.g., Moffet, 1968). However, it is convenient to assume a zero redundancy array in 
order to obtain an expression for AT which explicitly shows its dependence on the number of individual antennas in the 
synthesis array. Hence, let MN = n(n - l)/2 and, since n is likely to be large for any practical application in remote 
sensing from space, one may further simplify by writing MN = n 2 /2. Thus, in the case of a zero redundancy array of n 
antennas, Equations 29 can be written: 

T A 

distributed source 


T A 

sys syn 

n/B T n A„ 


(30a) 


AT =< 




T S) , Q c 

VB t n 


point source 


(30b) 
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Equation 30a is the expression commonly quoted in the literature on applications to radio astronomy (e.g., Napier et al., 
1983). Clearly, in either case (point source or distributed source), reducing the product, n A e , worsens the sensitivity. 
Hence, in applications to remote sensing from space there is a trade to be made between thinning the array to reduce the 
size and weight in orbit, and obtaining the best possible sensitivity. 

4) Hermitean Symmetry 

More data are available to create the image than the MN = n(n - l)/2 points discussed above. In particular, the 
signal < < x m x n > > at any baseline (m,n) is Hermitean: < < x m x n > > = < < x_ n x_ m > >* which can be seen 
from Equation 10 using the vector baseline L = md x x + nd > y. Hence, given measurements at MN baselines, the total 

A A 

number of points available for each dimension in the DFT are M = 2M + I and N =2N + 1 which are made up of the 
m = 1, 2, 3. . . M and n = 1,2,3. . . N actual measurements, the complex conjugate of these measurements and the 
measurement at zero baseline. Using all of these points in the DFT creates an image which is real (real values of bright- 
ness temperature). Using all of them results in a gain function G(q) which is real and symmetric. However, since the con- 
jugate points do not represent physically different measurements, the noise in these data is not independent. Hence, 
including them in the processor does not improve the sensitivity. These additional data points improve the resolution but 
not the sensitivity. 

IV. SELECTED CONFIGURATIONS 

Equation 26 and its special cases (Equations 29 and 30) provide the basis for the evaluation of several different vari- 
ations of aperture synthesis which have been proposed for remote sensing from space. Expressions for the sensitivity of 
the image will be derived in this section for the following configurations: a) a scanning real aperture antenna using a total 
power receiver (a special case of the correlation receiver); b) a representative synthetic aperture radiometer in which the 
antennas are arranged along the arms of “T;” c) a hybrid real/synthetic aperture sensor (called “ESTAR”) which has 
been proposed to make soil moisture measurements from space as part of NASA's Earth Observing System (Eos); and d) 
a variation in which individual antennas are arranged in a “cross” or “T” but instead of doing traditional aperture syn- 
thesis, the antennas are combined to form a linear array with multiple beams and the output from orthogonal beams are 
correlated to obtain resolution in two dimensions like a multiple-beam Mills cross (e.g., Mitman, 1988). Expressions for 
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the sensitivity of each variation will be obtained in this section and then the sensors will be compared for a hypothetical 
remote sensing application in the section to follow. 

A. Scanning, Total Power Radiometer 

The most common way a microwave radiometer is used to form an image is to scan the beam across the scene, 
spending t seconds at each spot. Systems with both electrical and mechanical scanning have been used for microwave re- 
mote sensing from space (e g., ESMR, SMMR, SSM/I). In such a system, each pixel is the result of a separate measure- 
ment and the sensitivity is just that of the receiver output. The receiver sensitivity can be obtained using the analysis 
in the preceding sections in the special case of zero baseline. Using the SNR as given in Equation 9, and following the 
same arguments which were used to obtain the sensitivity of the correlation receiver (Equations 10 through 17), one ob- 
tains the following expression for the sensitivity of the total power radiometer viewing an extended source: 

^ 

Vb T 

This is the conventional result for a total power receiver (e.g., Kraus, 1966) and as mentioned above, it applies to each 
pixel in the final image. 

Notice, that the sensitivity in this case (distributed source) does not depend on the antenna aperture. Also notice that 
Equation 31 differs from that obtained with aperture synthesis (Equation 30a) by the factor A syn / n A e which is the ratio 
of the effective area of the synthesized beam to the effective area of the antennas actually used in the array. Since n A c 
< A syn in any practical application, the synthesized beam will have a poorer sensitivity than a real antenna of equal effec- 
tive area if all other factors (i.e., B and r) are the same. However, as will be discussed below, the integration time, t, 
available per pixel can be very different for the real aperture and synthetic aperture sensors. In particular, because the 
synthetic aperture sensor doesn’t have to scan, the integration time can be much longer. 

B. Synthetic Aperture Array 

A synthetic aperture radiometer with a prescribed spatial resolution can be implemented with individual antennas in 
many different configurations. The only requirement is that the antenna positions provide the proper sampling in the spa- 
tial frequency domain. For example, the antennas could be arranged along the circumference of a circle, along the arms 
of a “cross,” along the arms of a “Y” or “T,” or even in 4 “random” array, and in each case obtain a synthesized 
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beam with roughly the same resolution. In designing an array for aperture synthesis, one has a great deal of flexibility to 
configure the system to accommodate the platform and to provide holes to avoid blocking the fields of view of other sen- 
sors. The sensitivity of a representative example will be discussed in this section to provide a reference. 

The example to be treated explicitly consists of identical antennas uniformly spaced along the arms of a “T.” This 
is a convenient example to analyze because the independent baselines which can be formed in such an arrangement map 
onto a Cartesian grid (x,y) where x is an integer multiple of the distance between antennas along one arm of the “T” 
(e.g., the top) and where y is an integer multiple of the distance between antennas along the other arm. If N x and N y are 
the number of antennas along each of the arms, then the number of independent baselines is the product N x N y , and the 
sensitivity obtained with this configuration is obtained from Equation 29a with MN = N x N y : 


T S ys 1 A syn 

AT = ■ 

y/rwr s/njt a 


(32) 


For the sake of simplicity, assume that each arm of the “T” is of equal length, L, and contains the same number of iden- 
tical antennas, each with dimensibn W (i.e., square patches W meters on a side). Then, to a first approximation A syn = 


L 2 and A e = W 2 and letting N x = N = N, the sensitivity becomes: 


T L 2 
AT = *y s 


vTlTF n w 2 

T I 

1 sys 

•JbT V2W 


(33a) 


(33b) 


where the last expression is obtained using the upper limit, N = L/W, to estimate the number of antennas. 

It may be possible to improve the sensitivity in some remote sensing applications by averaging pixels. In particular, 
the synthetic aperture radiometer forms an image of the entire field of view during each integration period (r seconds); 
and because the platform (e.g., a satellite in low earth orbit) is moving, the next image will be shifted slightly in the 
direction of motion. If the integration time is chosen so that this shift is just one pixel, then each pixel will appear in 
several successive images (N times in the example being discussed here). Each view will not be at the same incidence an- 
gle, but if this difference will not degrade the data, then it is possible to use these multiple looks to further reduce the 
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noise. Assuming that each look is independent, one obtains an improvement of 1/ VN", and using N = LAV the sensitivity 
(Equation 33) becomes: 

T f ~ 

sys 

AT = —i / (34) 

VB t V 2 W 

At first glance, the sensitivity of this synthetic aperture radiometer, even in the optimum form given in Equation 34, 
is poorer than is obtained with a scanning, total power radiometer (Equation 31); however, this is misleading because the 
integration time available to form an image in the two cases is very different. In particular, in order for the real aperture 
radiometer to form a map comparable to that obtained with the synthetic aperture radiometer (i.e., a map consisting of N 
x N contiguous pixels), it must complete a scan in t'/N seconds where t' is the time the spacecraft is over one pixel. The 
integration time available for aperture synthesis (Equation 34) is t' and the time available for the scanning, total power 
radiometer (Equation 31) is t'/N. Using t' = 1/N = W/L in Equation 31 yields a result which differs from Equation 34 
by only 1/V2 . Hence, in this case, synthetic and real aperture sensors with comparable system temperatures and band- 
widths provide comparable radiometric sensitivity. 

C. Hybrid Real and Synthetic Aperture Array (ESTAR) 

The sensitivity obtained using aperture synthesis (Equation 30a) depends inversely on the actual collecting area (nA € ) 
employed in the array and also inversely on the time-bandwidth product of the measurement. In remote sensing applica- 
tions in space, the integration time is limited by the platform motion (about 7 km/sec in low earth orbit) and the band- 
width is limited by RFI and requirements on the field of view of the instrument (Thompson and D’Addario, 1982). 

Hence, in general, a trade has to be made between the sensitivity which must be achieved for a given measurement to be 
useful and the reduction in antenna hardware (i.e., nA e ) which can be obtained. ESTAR is one such compromise which 
has been proposed to make measurements of surface soil moisture practical from space (Murphy et al., 1987; Le Vine et 
al., 1989). The idea is to increase sensitivity by doing aperture synthesis in only one dimension. ESTAR achieves spatial 
resolution in the along-track dimension by employing long stick antennas (real apertures) aligned parallel to the motion of 
the spacecraft and achieves resolution in the other (cross-track) dimension by means of aperture synthesis using pairs of 
the stick antennas. Using aperture synthesis substantially reduces the number of stick antennas required (Le Vine et al., 
1989; Swift et al., 1986). 

The sensitivity of this configuration is obtained directly from Equation 30a where A e is the area of the stick anten- 
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nas and n is the number of sticks employed. Letting the sticks be L meters long and W meters wide, and using the ex- 
pressions A e = WL and A syn = L 2 to approximate A e and A syn , one obtains: 


T L 

sys 

vETT n W 


( 35 ) 


This is the sensitivity of each pixel in the image formed by this hybrid sensor. The ratio, L/(nW) , is the “fill factor’ ’ of 
the array: it is the fraction of the total available area which is occupied with receiving antennas. The goal is to reduce the 
fill factor as much as possible to save weight while at the same time achieving the required sensitivity. If one assumes 
that r is fixed by the velocity of the platform and B is limited by fringe washing, then the ratio L/nW determines the 
amount of thinning that can be achieved and still meet the sensitivity requirement of the measurement. In the application 
to soil moisture, this array could be more than 80 percent empty space (e.g., Le Vine et at., 1989). 


D. Correlation Cross Antenna 

Another sensor system to which the formulas above apply is the variation of the Mills cross described recently by 
Milman (1988). In this system, two linear arrays are arranged orthogonally in the form of a “cross” or “T.” 

The antenna elements in each arm are used to form a collection of contiguous beams with resolution in the dimension per- 
pendicular to that arm. The intersection from a pair of beams (one from each arm) determines the resolution 4 ‘cell” and 
the signal from this cell is detected by correlating the output from this pair of beams. For example, let X t be the i-th 
beam formed with an array along the x-axis and Y J be the j-th beam formed with an array along the other arm. Then, in 
the proposed system, the signal from the intersection of these two beams is detected by forming the product 
<X i (t)Y j *(t)> in the processor. Radiation from the entire scene can be mapped using contiguous beams and forming the 
product of all pairs (i,j). 

Although this system could be implemented with a collection of beam forming networks and correlators (Milman, 
1988), it can also be implemented as a special case of aperture synthesis in which the individual antennas in the synthesis 
array are aligned along the arms of the cross. To see this, first note that the correlator required to form the product is the 
same as indicated in Figure 1 (correlation receiver) with the two arms (Xj) representing the signal from two beams. Con- 
sequently, the result of multiplying the two beams together is the same as given in Equation 2a, but with S { j(v) represent- 
ing the spectrum formed from the product of two beams. It is the Fourjer transform of the correlation function 
<Xj(t)Yj(t + r). Now, to compare with aperture synthesis, assume that each arm consists of 2N + 1 identical, uniformly 
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spaced elements whose output voltage (in the frequency domain) is V n (v) where n = 0,±1,±2,±3 ... ±N. Then the i-th 
beam is formed by adding an appropriate phase to each element in the array (e.g., Jordan and Balmain, 1968). In the fre- 
quency domain this can be written: 

X,( v ) = output of the i-th beam (36) 

N 

= V x„( v ) ex P(i n k d qJ 

n = -N 


where k = 2 tt v/c, d is the spacing between elements and q xi is the direction cosine of the vector in the direction of the 
main lobe (main beam) of the i-th beam measured relative to the axis of the array ( - 1 < q xl < +1 and q xi = 0 at 
broadside). Now, the correlation function for two beams can be written in terms of its frequency domain representation as 

follows: oo 

<X,(t + r)Yj*(t)> = [ <X t (v)Y*(v)> exp(-j 27 rvT) dv (37a) 


and using Equations 36 and the assumption that < V xn (v)V ym *(v')> = < V xn (v)V ym *(v) > <$(v - v') and re-ordering 
the sum on n, one has: 

N N 

<Xj(v)Yj*(v)> = < V xn( v ) V yn,*( v )> exp[ -jkd(nq xi + mqj] (37b) 


= - N n = - N 


The average <x i (v)X J (v)> is the spectrum Sy(v) in Equation 2a when the inputs to the correlator (Figure 1) are the 
beams formed by the two arms of the cross. Notice, that it has the form of a discrete Fourier transform of the average 
< V xn( v )V ym (v) > , which in turn, is the spectrum Sy(v) in Equation 2a when the inputs to the correlation receiver are the 
signals from the individual antennas in each arm. Hence, Equation 37b represents the signal processor required to form an 
image when employing aperture synthesis. Thus, the process of forming and correlating a pair of beams is just an alterna- 
tive way of implementing the “processor” in Figure 1 required in aperture synthesis. As long as the two arms of the 
cross intersect (i.e. have an element in common) then the sum in Equation 37b involves all possible pairs of antennas and 
the image formed is identical to that obtained with aperture synthesis when employing the same configuration of elemental 
antennas. 
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It follows that Equations 16 and 26 for the sensitivity of a synthetic aperture imaging system also apply to this con- 


figuration. In the case when the scene is extended. Equation 26 reduces to Equation 30a and one has: 


AT = 


sys 


syn 




(38) 


/ M N A. 


This applies to the cross configuration being considered here where T sys , B and t are the system noise temperature, band- 
width and integration time of the correlator used for each pair of beams; A c is the collecting area of each elemental anten- 
na; A syn is the effective area of the beam formed by each arm of the cross; and MN is the number of independent pairs 
of individual antennas in the cross. 

In order to facilitate comparison with the ESTAR configuration described above, assume arms of length, L, and 
width, W, and suppose that there are N = L/W antennas in each arm. Then, A syn = L 2 and A e = W 2 , and there are 
N 2 /2 independent antenna pairs in this array. (Only the baselines formed between one arm and half of the other arm of a 
cross, “+”, are independent.) Thus, with MN = N 2 /2 = L 2 /(2W 2 ) Equation 38 can be written: 


T s s L 

AT = -== — (39a) 

VbT W 


However, it may be possible to obtain better sensitivity by averaging the pixels in successive images as done above (Sec- 
tion IVB). This can be done on a moving platform by adjusting the integration time so that the shift due to platform mo- 
tion during the integration is just one pixel. Then each pixel will be appear in several successive images (N times). 
Assuming that each look is independent and that the difference in incidence angle will not degrade the data, one obtains 
an improvement of 1 /Vn = VW/L and Equation 35 becomes: 

T s s L 

AT = -=. ~ (39b) 

VBt W 


In summary, it should be clear from this discussion that one could implement this system either with a collection of 
beam forming networks and correlation receivers, or by forming all the n correlation pairs first and then doing the beam 
forming (a DFT). In the latter case the system is clearly a conventional form of aperture synthesis. 


V. COMPARISON 

In order to compare the performance of the several configurations described above, consider as a specific example 
the problem of measuring soil moisture from space. Imagine a sensor in polar orbit as part of the proposed Earth 
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Observing System” (Butler et al., 1988). Preliminary studies call for an instrument at L-band (1.4 GHz) with a spatial 
resolution of about 10 km (Murphy et al., 1987; Le Vine et al., 1989). At the orbit proposed for the Eos (824 km) this 
means an antenna aperture on the order of 20 m on a side, clearly a candidate for aperture synthesis in some form. It will be 
assumed that the platform is moving at 7 km/s and that a swath of 1000 km (global coverage every 3 days) is required. 
The integration time will be chosen to be the time available per resolution cell. That is r = 10 km - 5 - 7 km/sec for the 
synthetic aperture sensor and r = (10 km -5- 7 km/sec) + 100 for the scanning, real aperture radiometer (where 100 is 
the number of resolution cells per scan line). A bandwidth of 10 MHz will be assumed for the synthesis arrays (a 
reasonable value for which fringe washing should not be a problem and digital processing can be readily implemented) 
and a bandwidth of 30 MHz has been assumed for the total power radiometer (this is the full bandwidth available at the 
1.4 GHz radio astronomy band and, although a larger bandwidth is technically feasible, a system with a larger bandwidth 
is likely to encounter severe RFI problems at L-band). A system noise temperature T sys = 500K will be used for all sen- 
sors. This is a conservative estimate, and lower system temperatures likely can be achieved. 

Table I is a comparison of the sensitivity of each of the candidate configurations discussed above. The first row in 

the table gives the general theoretical expression for the sensitivity in the image using the formulas derived above for the 
case of a distributed target. The second row gives the sensitivity for the specific configuration under consideration. The 
synthetic aperture configurations being considered are: 1) a representative configuration with individual antenna elements 
uniformly spaced along the arms of a “T”; 2) the “cross” configuration; and 3) the hybrid synthetic and real aperture sensor, 
“ESTAR.” Expressions for the sensitivity of the former two configurations are given both with and without additional 
pixel averaging. The total power radiometer is assumed to have an antenna 20 meters on a side and the length of the ar- 
ray in each of the synthetic aperture configurations will be assumed also to be 20 m to achieve approximately the same 
resolution with all of the sensor configurations. The width of the array has been assumed to be about one-half a 
wavelength (0.1m) and it is assumed that each elemental antenna in the array is W =0.1 m on a side. This permits a 
minimum spacing between antennas of one-half a wavelength in all the synthesis arrays. An improvement in sensitivity 
could be achieved by using bigger elemental antennas (larger W); however, when W > X (approximately) the grating 
lobes are no longer in imaginary space. The number of antennas in each array is L/W except for ESTAR. In the ESTAR 
configuration, aperture synthesis is done only in one dimension, and in this configuration the number of sticks has been 
chosen to be the number needed in a minimum redundancy array to achieve spacings which are an integer multiple of 
X/2. In practice, more antennas would be needed, but no additional independent spacings would be achieved. The cross is 
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assumed to be in the form of “ + Thus, one only uses N/2 elements in one of the arms which explains the factor of 2 
in the formulas appearing in the Table. The “ + ” will also affect resolution somewhat, but this is being ignored. The 
cross could be implemented in the form of a “T** in which case, the formulas are exactly the same as for the synthetic 
aperture array having the same form. 

The third row in Table I gives a numerical value for the sensitivity obtained using the parameters shown at the bot- 
tom of the Table. Notice, in particular, that when pixel averages are employed in the ‘‘cross” and “T" configurations, 
the sensitivities of both the real aperture and synthetic aperture are about the same. However, without pixel averaging, 
only the “ESTAR” configuration offers a sensitivity, AT, close to that achieved with a real aperture, scanning radi- 
ometer. 

VI. CONCLUSIONS 

General expressions have been derived for the sensitivity of a radiometer which uses aperture synthesis to form an 
image. Among the interesting conclusions is that the sensitivity of the individual correlation receiver, of which such a 
radiometer is comprised, depends on both the baseline (distance between antennas) and on the nature of the scene (e.g., 
point source compared to distributed source). The same is true of the image formed with aperture synthesis from a collec- 
tion of such receivers. The sensitivity can be different for different pixels in the image and in addition, can depend on the 
nature of the source. This result is different than in the case of a conventional imaging radiometer system (i.e., a scan- 
ning total power radiometer) where a single receiver and a single antenna are used to form an image. In this case, the sen- 
sitivity can be made independent of the antenna and scene by an appropriate definition of the system noise temperature 
(e.g., see Equation 18). 

In the situations likely to be encountered in microwave remote sensing of the earth from space, the expressions for 
the sensitivity given above can be simplified (e.g., assuming a rtiinimum redundancy array with a large number of ele- 
ments and a distributed source). In this case, the sensitivity obtained with aperture synthesis is proportional to that 
obtained with a total power radiometer of the same system temperature, bandwidth and integration time (i.e., T sys /\/T5r). 
The proportionality constant is the “fill" factor, A syn /n A tf , which is the ratio of the effective area of the synthesized 
antenna to the actual collecting area employed in the array. 

The advantage of aperture synthesis is that it can achieve spatial resolutions equivalent to total power radiometers 
with large effective collecting areas (i.e., equal to A syn ), but using small antennas (A c ). The reduction in sensitivity that 
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this entails can be restored because the synthetic aperture system does not need to scan and it collects energy from many 
independent antenna pairs. The comparisons given in Table I indicate that very comparable sensitivity can be obtained 
with several synthetic aperture systems which have very substantial reductions in collecting area. Although the comparison 
applies specifically to the special case of soil moisture measurements from space, similar arguments should apply to other 
remote sensing applications as well. 

There are many possible ways to implement aperture synthesis; and each configuration must be evaluated for its sen- 
sitivity. General formulas have been given in this papier which can apply to a wide variety of situations and antenna con- 
figurations. 
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Figure 1. Schematic illustrating the steps involved in aperture synthesis. 
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Table I. The sensitivity obtained with several types of aperture synthesis and with a total power radiometer. 
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APPENDIX A 


THE MEAN AND VARIANCE OF THE OUTPUT OF A 
CORRELATION RECEIVER 


In this appendix the mean and variance of the output y(r) of the correlation receiver shown in Figure 1 will be com- 
puted. In doing so, the following notation will be used: 


Sij(v) 


oo 


1 


Ry(t) exp(j2irvt) dt 


- OO 


(Al) 


where R^t) = <r i (t)r J (t + r)> is the correlation function of the input voltages, r(t), and the S^v) are the power spec- 
tra associated with these correlation functions. It will be assumed in all that follows that the r j(t) are stationary random 
processes. In describing the Filters, lowercase symbols [e.g., h(t)] will be used to denote functions of time, and capital 
letters [e.g., H(v)] will denote their Fourier transform. The symbol will be used to denote a convolution, except 
when used as a superscript when it will denote the complex conjugate. 


A. MEAN 

Referring to Figure 1 and letting <y r (r)> be the mean output of the radiometer when the input is r, one obtains: 


<y,j( r )> = <e,ej> * h 


LP 


(A2) 


= < [r, * h,] [r j * hj] > * h LP 


oo 

-i 


R 1J (t ' - t") hid'” - f) hj(t " ' - t ") h LP (t - t" ') dt 'dt "dt " 
Now, replacing all the terms in the integrand by their Fourier transform, one obtains: 


oo 

<y jj (r)> = H lp ( 0) | S^v) H^v) H^fv) dv 


(A3) 
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B. Variance 


The variance of the receiver output is: Var[y] = <y 2 > - <y>. Since the mean value, <y>, is given above 
(Equation A3), it is necessary only to compute <y 2 >: 

<yij 2 (r)> = < [z,j * h L p] [Zjj * h Lpl > (A4) 

00 

= j S zz (v) H lp 2 (v) dv 
- 00 

where S zz (v) is the Fourier transform of the correlation function R zz (t) = < Zy(t)z t j(t + t)>; and the second line is ob- 
tained by substituting for the filter in terms of its Fourier transform as done above (e g.. Equations A2 and A3). Sub- 
stituting for Zyft) in terms of the signals in each arm of the receiver (i.e., Zy(t) = e | (t)e j (l); Figure 1), one obtains: 

R zz (r) = < [e | (t)e j (t)] [e,(t + rfe/t + r)J > (A5) 

= <e,(t)e,(t + r)> <e j (t)e j (t + t)> 

+ <e i (t)e j (t)> < e j (t + r) ej (t + t)> 

+ <e,(t + r)e j (t)> <e i (t)e j (t + r)> 

where the last step is obtained by assuming that the e,(t) are Gaussian random processes. Now, following the same 
procedure as used above to calculate the mean, one obtains the following result for each of the pairs above: 

< e^Oe^t + t)> = < [r,(t) * h ( (t)] [r^t + r) * hj(t + r)] > (A6) 

oc 

= | Sy(v) Hj(v) Hj*(v) exp(-j2iTVT) dv 

<x> 

Equation A6 applies to any combination (ij): that is, for <e i ej> as well as < e^j > and < e j e j > , etc. Now, taking the 
Fourier transform of Equation A5 using the results of Equation A6 above, one obtains: 


oc 

SJv) = 6(v - 0) { | Sy(v) H,(v) Hj*(v) dv } 2 (A7) 

— ao 

+ [ S jj (v)Hj(v)H i *(v) ] * [ Sy(v)H j (v)H j *(v) J 
+ [ S 1J (v)H l (v)H J *(v) ] * [ Sy(v)Hj(v)Hj*(v) ] 

= S“Jv) + S b zz (v) + S'Jv) 
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Substituting this result (Equation A7) into Equation A4 and comparing with Equation A3, one obtains: 


<y,j 2 (0> = <yij(r) > 2 + l [S a zz (v) + S b zz (v)] H 2 LP (v)dv 


It follows from Equation A8 that: 

Var [yy(r)] 


oo 



[S b zz (v) + S c zz (v)] H 2 Lp (v) dv 


= [ S b zz (0) + S c zz (0) ] 


OO 



H 2 LP (v) dv 


(A8) 


(A9) 


The last step in Equation A9 follows because it has been assumed that H LP (v) only passes signals in a narrow band very 
close to v = 0. Finally, substituting for the S zz (0) from Equation A7, one obtains: 


Var[y (J (r)] 


00 


! 


— oo 


[S 2 ,j(v) + S jj (v)Sjj(v)] H 2 i (v) H 2 j(v) dv H 2 lp (v) dv 


(A 10) 


The mean and variance of the output of the correlation receiver shown in Figure 1 are given in Equations A3 and 
A 10 for a reasonably general set of assumptions. They simplify somewhat in the important, special case when the system 
transfer functions H|(v) are very narrow compared to the spectra of the input signal, S l} (v), which is the situation encoun- 
tered in systems designed for microwave remote sensing of the earth from space. In this narrow bandwidth case, the sig- 
nal spectra can be factored out of the integrals and one obtains: 


< y 


oo 

,(r) > = Sy(v,) H lp ( 0) j H,(v) H *(v) 


dv 


(All) 


OO OC 

Var [ yjj (r)] = [ S 2 ij (v„) + S^v^S^v,,) ] . J H 2 i (v) H J /v) dv | H 2 LP (v) 


dv 


(A 12) 


where v 0 is the center frequency of the Filters, H } (v). Equations A 1 1 and A 12 are the ones used in the text to compute the 
signal-to-noise ratio of the receiver output. 




APPENDIX B 


POWER SPECTRUM OF THE ANTENNA OUTPUT VOLTAGES 

It is the objective of this appendix to derive an expression for the spectrum of the antenna output voltages (Figure 1) 
which are the inputs to the arms of the correlation receiver. In the notation employed in the text this spectrum has been 
denoted, < <XjXj> > , where the double brackets < < > > denote a Fourier transform and X l X j is a shorthand notation 
for the correlation function < + r)> . (In the text, a distinction was made between the response, X°, to a uniform 

scene and the response, x i? to fluctuations about this background. In this case one has: = X° 4- This distinction is 

not being made in this appendix and Xj is the response of the antenna to an arbitrary scene.) Referring to Figure 1, the X t j 
are the same as the voltages, at the antenna output terminals; and the discussion in this appendix will be in terms of 
the voltages, V(r if t). 

In this appendix, let V(rj t) denote the time-dependent voltage at the output terminal of the antenna at position r^ and 
using the same notation as adopted in Appendix A, let S^(v) be the Fourier transform of the correlation function Ry(r) = 

< \(r i (t + r))V*(rj t )> . In the special case when the voltages are stationary random processes, one has S^v) = 

t 

<V(r i ,v)V*(rj,v)> where the V(rj,v) are the Fourier-Stieltjes transforms of the antenna voltages. 

In order to obtain expressions for the V(r|,v), imagine two parallel planes, the scene plane which is the source of 
the radiation (e.g., the surface of the earth) and a second plane called the image plane at a height H 0 above the scene 
plane. The image plane contains all the interferometer antenna pairs. Let the coordinate origin be in the scene plane and 
let r' = + x'x + y'y + z'£ denote a position in the scene plane (z'=0) and r- — x 4 x + y-y + ZjZ denote the position of 
an antenna in the image plane (z l = H 0 ). The geometry is illustrated in Figure Bl. 

It is desired to obtain an expression for the output voltage V^v) from the antenna at jj in terms of the electric 
field, E s (r/,v), on the scene plane. To begin, consider radiation from a small patch Ax' Ay' at i' in the scene plane. The 
electric field AE(rj/r \v) radiated from this small patch to the antenna at r { can be obtained from the vector Helmholtz 
equation (by assuming that the patch is a small aperture in an otherwise opaque screen on which the tangential compo- 
nents of the fields are zero). Following standard procedures (e.g., Tai, 1971), one obtains: 

AEfc/r'.v) = - { [n x V x E s (r»] •G(t , /t) + [nx E s (r',v)]. V x G(r'/r.) } Ax' Ay' (Bl) 

where n = z is the unit vector normal to the scene plane and G(r i /£ f ) is a dyadic Green’s function satisfying the free 
space wave equation in the region z'> 0. A convenient choice for GijiJl!) is (Tai, 1971; Section 18): 
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G(r i /r') = [1 - 1/k 2 VV ] [ g 0 (r i /r') - g 0 (r i /r') ] - 2zz ^(r/r 0 (B2) 

where I is the unit dyadic, g 0 (r/r/) = exp(jk | r-r' | )/4ir | r-r' | is the free space scalar Green’s function and L = L~ 
2z£is the “image” of a point at r behind the scene plane. The following properties apply to this Green’s function when 


the source point £ ' is on the scene plane: 

n x G(r 7r) = 0 (B3a) 

V' x Gd'/r) = V' x I g 0 (r'/r) - V' x\ g 0 (r'/r) (B3b) 

(n x E).V'xG(r'/r) = 2 (jk - 1/R) gJ&L’) [ V R X (g x z) ] (B3c) 

| - 2 z z (B3d) 

Now, noting that [n x (V x E) ]• G = -(V x E)*[n x G] and using the preceeding properties of the dyadic 
Green’s function, one may write: 

ABd/r'.v) = 2 (jk - 1/R) g 0 (r | /l') [ V R x (E s x z) ] (B4) 


The expression in Equation B4. is the signal radiated from the patch on the scene plane to an antenna at r . The 
response of an antenna to radiation from the small patch can be expressed in terms of the voltage transfer function of the 
antenna, A(r/r',v), as follows: 

AV(r/r',v) = AE(r/r',v) • A(r/r',v) (B5a) 

where the voltage transfer function is a vector quantity to account for the polarization of the antenna (Collin and Zucker, 
1969) and can be obtained from the reciprocity theorem for antennas. One can show that A(r/l',v) • A*(r/r \v) = 

A e P n (r/r 'J where A e is the effective area of the receiving antenna and P n (r/r 0 is the normalized power pattern of the an- 
tenna (Kraus, 1966). 

To obtain the response of the antenna to all sources, it is now necessary to sum over all of the scene plane. In the 
limit of infinitesimally small patches, the sum approaches an integral and one obtains: 

V(lj,v) = ^ AEfr/r'.v) • A(i/r',v) dr' (B5b) 

= 2 j E 0 (r',v).A(r i /r') [jk - 1/RJ [z,/R, ]g 0 (r/l') dr' (B5c) 

A A 

where the radiated electric field in Equation B5c has been expressed in terms of it components, E 0 (r',v) = e h h + e v v, 
along horizontally and vertically polarized unit vectors (fi, v). These unit vectors are defined locally (i.e. at each point r') 
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with respect to the plane formed by the normal to the surface (z) and the vector, R t = r { - v\ from the antenna at r i to 
the source point at r\ The unit vectors are: 

h = (z x V R.)/sin(^|) (B6a) 

0 = V R t x (z x V R i )/sin(fl i ) (B6b) 

where sin(0 t ) = | z x VRj | is the angle between the line-of-sight and the z-axis. The projections of the electric field 
onto these unit vectors (e h , e v ) can be expressed in terms of the field on the surface, E s (r\v), as follows: 

^ = E s *h cos(d { ) (B6c) 

e v = [ (E s *v) - (E 1 .VR i )tan(ff i ) ] cos(0j) (B6d) 

In Equation B4 the factor, cos (0 { ), appearing above has been factored out and appears explicitly in the form cos(^)= z l /R 1 
where Rj is the magnitude of the vector R t and z { = H 0 (Figure Bl). Hence, in Equation B5b one has e h = E s *fi and e v 
= £ s .v - (E^VR^tantf,). 

Now, using these results and assuming that the distance between the image plane and scene plane is much larger 
than the wavelength of any radiation in the passband of interest (i.e., kR t > > 1 for all k of interest) and that the fields 
on the surface are (spatially) uncorrelated, one obtains: 

Sy( v ) = < V(r,,v) V'dj.v) > = (k 0 /2ir)- | (H^R^) 2 < A(r i /r',v).T(r',v).A*(r J /r,v)> p exp(jk[R, - R^) dr’ (B7) 

where the dyadic intensity. I(r\v), has been introduced to account for the polarization of the incident radiation. It is 
defined as follows: 

7(r\v) = 

where I ah = < e a (r',v) e b *(r',v) > . The brackets < > in the definition for I ah denote a statistical average over the spa- 
tial variable of the random process, and the brackets < > p in Equation B7 indicate an average over polarization states 
(Collin and Zucker, 1969). 

In the practical cases which are of interest in remote sensing from space, the distance between the scene and image 
planes is much larger than the distance between the receiving antennas. In this case one may approximate R^ in the 


IhK h h 


Ik. V h 


V h v 


I A A 

I., V V 


(B8) 
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denominator in Equation B7 by RjRj = R 0 2 , and for the factor R, - Rj appearing in the exponential, one may write 
(Appendix C): 

R, - Rj = | Ij - r'| - 1 Ij - ii | = (x, - xp cos(a x ) + (y, - y^ cos(a y ) (B9) 

where cosfcx ) are the direction cosines of the vector from the source point at £_ in the scene plane to the origin in the 
image plane. Using this approximation, one obtains: 

Sy(v) = f [k 0 H 0 /2 R 0 2 ] 2 < A(Tj/r',v)» T(r',v)«A*(rj/r',v) > p exp[j2jr(i7 x cosa x + r? y cosa y )] dr' (BIO) 


where 


x, - Xj 


% = 


y> - yj 


\ = 


(Bll) 


R o = |H 0 S -I'l 

It is convenient to do the integration in Equation BIO in terms of the direction cosines cos(o: x y ) rather than in terms 
of the coordinates r Thus, making the following change of variables: 

q x = cos(a x ) = x'/R^x'.y') 
q y = cos(a y ) = y7R 0 (x',y') 

and using the Jacobian of this transformation dq x dq y = [H 0 /R 0 2 ] 2 dx 'dy one obtains: 


(B12) 


+ 1 


S (J (V) = (l/X 2 ) 


r* 

< a.t.a' 

t 


A*> exp(j2ir2*g) d£ 


(B13) 


-l 


where 3 = ij x x + ij y y and ^ = q x x + q y y . 

A special case of importance is that in which the antennas are identical and the radiation is unpolarized. Ignoring 
the slight difference in location of the antennas (i.e., assuming that Afr/r'.v) = Afrj/r'.v) ), one obtains: 


< A-I-A*> p = (1/2) A e P„(q) yq) 


(B14) 


where P n (£|) is the normalized power pattern of the antenna in the direction (j (Kraus, 1966), A e is its effective receiving 
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area, and I^q) is the spectrum of the electric field on the scene plane with polarization p: I pp (£) = < [p*E 0 (r,v) ] 2 > p . 
Hence, in this case: 

+ 1 

Sjj(v) = (A e /2X 2 ) j P n (q) I pp (q) exp(j2in,*q) dq (B15) 

In Appendix E it is shown that I pp (g)/(Z 0 \ 2 ) is the brightness in the direction Hence, 

+ i 

S 1J (v)/Z 0 = (A e /2) j P n (£) B(£) exp<j27r r a) d£ (B16) 

-I 

where the definition B(£) dq = B(ft) dQ has been made (dft = sin(0)d0d<£ denotes the “solid angle” in spherical coor- 
dinates). It is also shown in Appendix E that B(q) = B(cj) / Vl - q x 2 - q y 2 . Hence, one may write: 

V p„(a) B (a) 

S ;: (v)/Z„ = (A^/2) l exp(j27ri7»q) dq (B17) 

J VI - q x 2 - q y 2 
- ! 

where B(q) is now the conventional definition of brightness and Equation B17 is the same as has appeared in the literature 
(e.g., Napier et al., 1983). Also, notice that S iJ (v)/Z 0 is the fraction of the power incident on the receiving antenna (per 
unit bandwidth) which is available at the antenna output terminals. It is the power which can be delivered to a load 
properly matched to the antenna (e.g., Collin and Zucker, 1969). 

REFERENCES: 
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APPENDIX C 


THE SAGITTAL APPROXIMATION 


As defined in the previous discussion (Appendix B), r t — x, x 4- y, y +H C z is the position vector of an antenna in 
the image plane, r' = x’x + y’y is the position vector of an arbitrary point on the scene plane, and Rfr/lO = I i — r is 
the vector between the point r’ on the scene plane and the antenna. In this appendix an approximate expression for 

R(r/r ^ - R(lj/r ^ will be derived where R = | R | . 

Letting R 0 = H 0 z-r' and assuming that r ti • r (J < R 0 2 where r tj = x,x + y,y, one can write: 


R(r/i') - Rd/iO = I l - I' I - lij - 1' 


= I Bo + in I 


.Bo -IJ 


= VR 0 2 + r tj 2 - 2 R 0 .r ti - VR 0 2 + r 2 _ 2 R o ._ r 


(Cl) 

(C2) 

(C3) 

(C4) 

(C5) 

(C6) 

(C7) 


= R 0 Vl + (r ti 2 /R 0 2 ) - 2 R 0 *r ti /R 0 2 - R 0 Vl + (r tJ 2 /R 0 2 ) - 2 R 0 *r tJ /R 0 2 
= (r ti 2 - r tJ 2 )/2R 0 + R 0 .(f tj - r ti )/R 0 
= R 0 -(r,j - I,,) /R o 

= (x, - Xj) cos(a x ) + (y, - y^ cos(a y ) 
where cos(a ) are the direction cosines of R 0 (a vector from the source point to the origin in the image plane) with 
respect to the coordinate axes in the image plane [cos 2 (a x ) + cos 2 (a y ) + cos 2 (a z ) = 1 ]• 

In order to use this approximation in the exponential in Equation B7, the error made in replacing Equation C4 by 
Equation C6 must be small. Hence, one must have k(r ; 2 — rj 2 )/R 0 < <TT. As an upper limit on this error one can replace 
( r 2 _ v 2 ) w j,h l 2 where L is the maximum spacing between antennas. Hence, the approximation requires that the size of 
the array in the image plane (squared) divided by the minimum distance to the scene plane be much less than a 
wavelength. This is a condition easily met in microwave remote sensing of the earth from space (e.g., even in an extreme 
L is not likely to be more than a few 10 ’s of meters, but even in a very low orbit, R t , is on the order of 300 km or 


case 


more). 
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APPENDIX D 


COORDINATE TRANSFORMATIONS 


In this appendix the transformation between the Cartesian coordinates (x', y\ z 1 ) and the direction cosines (q x — 
cos(a x ); q y = cos(a y ); q z = cos(a z )) is established. Letting R(x',y',z') = | r' | , one has: 

q x = x' / R 

q y = y'/R (Dl) 

q z = z' / R 

and q x 2 + q y 2 + q z 2 = 1. In the case at hand, z' = H (> is a constant and of interest is the relationship between dx'dy' 
and dq x dq y . Using the Jacobian of the transformation to establish the required relationship, one obtains: 


dx' dy' 

dq x dq y = dx ' dy ' (D2) 

dq y dq y 

dx ' dy ' 

= (H 0 / R 0 2 ) 2 dx' dy' (D3) 

In establishing Equation D3 the following relationships have been used: 

dq x /dx = [1 - (x/R) 2 ]/R (D4a) 

dq x /dy = - (xy)/R 3 (D4b) 

dq y /dx = — (xy)/R 3 (D4c) 

dq y /dy = [1 - (y/R) 2 ]/R (D4d) 


The relationship between the direction cosines (q x ,q y ) and the spherical coordinates (6,0) has also been used in the 
text. The coordinate transformation is: 


q x = cos(a x ) = sin (0) cos(^) 
q y = cos(a y ) = sin (0) sin (</>) 


(D5) 
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and the Jacobian of the transformation is: 

‘K dq x 

d0 d<t> 

Jacobian = = sin(0) cos (<t>) 

d Q, dq y 

d 0 dc/T 

Thus, one has 

dq x dq y = sin(0) cos(0) d0 d</> 

— cos(0) dO 

where dQ = sin(0)d0d<£ is the differential “solid angle.” 


The inverse transformation is: 



(D6) 


(D7) 


(D8) 
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APPENDIX E 


RELATION SHIP BETWEEN BRIGHTNESS AND ELECTRIC FIELD INTENSITY 


In the notation commonly used in radio astronomy (e.g., Kraus, 1966), the radiation from a source is defined in 
terms of the source “brightness.” The brightness is defined to be the power per unit area arriving at the receiver per 
steradian from the direction (0,0) and in a narrow frequency interval, Av. In particular, B(0,0)sin(0)d0d0 = B(Q)dD is 
the power per unit area per unit of bandwidth arriving at the receiver. The radiation received from a thermal source can 
be treated as if it were produced by electromagnetic fields on the surface of the source just as radiation from an antenna 
can be treated in terms of the equivalent electromagnetic fields on the aperture of the antenna. The objective of this 
appendix is to relate B(I2) to these equivalent electromagnetic fields on the surface of the source. 

In an electromagnetic system, power is obtained from the Poynting vector, P — E x H. When the observer is very 
far from the source, the electromagnetic fields are nearly plane waves and the magnitude of the Poynting vector is ap- 
proximately: P = E-E/Z 0 where Z () = V*i 0 /€ 0 is the characteristic impedance of free space. Assuming that the fields are 
ergodic, stationary random processes (in the time coordinate), one obtains: 

< P(r,t) > = < E(r,t)*E(r,t) > / Z G 

00 

= ^ P(r,v) dv (El) 

— oo 

where P(r,v) is the Fourier transform of the autocorrelation function < E(r,t)*E(r,t-br ) > evaluated at the position r It 
can also be expressed in terms of the Fourier-Stieltjes transform, E(r,v), in the form: P(r,v) — <E(r,v)*E*(r,v)>/Z 0 . It 
follows from Equation El that the time average power in a frequency band, Av, which is narrow compared to PQ;,v) is 
approximately P(r,v) Av (e.g., Papoulis, 1965). 

The spectrum P(r,v) of the radiation at the receiver at £ due to a patch of surface at £' can be calculated using the 
procedure outlined in Appendix B (letting i = j). Assuming that the radiation separates in space and time coordinates and 
is spatially uncorrelated, one obtains: 

P(r,v) dr' = [H n /R 2 (r/r , )l 2 [l/X 2 ] <E t a\v)-E,(r\v)/Z 0 > p dr' (E2) 
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Comparing this expression with the definition of brightness, one obtains: 


B(fl) dU = P(r,v) dr' = [1/X 2 ] <E t (r',v)*E l (r»>/Z 0 > p [H 0 /R 2 (r/r')] 2 dr' (E3) 

Changing from Cartesian to spherical coordinates on the right-hand side in Equation E3, one has [H 0 /R 2 (r/r')] 2 d£' = 
cos(0) dft. Hence, the relationship between the brightness and the electromagnetic fields on the surface is: 

B(Q) = [1/X 2 ] [<E t (fi,v)*E t (Q,v) > p /ZJ cos(0) (E4) 

The brightness can also be expressed in terms of the direction cosines (q x , q y ). Defining an alternative brightness, B(^), 
and noting that (Appendix D) dq x dq y = (H 0 /R 0 2 ) 2 dx'dy', one obtains: 

B(q) dq = P(r,v) dr' = [1/X 2 ] [ < E t (q,v).E t (q,v) > p /Z 0 ] dg (E5) 

Finally, note that 

B(q) = B(fi)/cos(0) = B(q)/Vl - q x 2 - q y 2 (E6) 


REFERENCES: 

Kraus, J.D. (1966), Radio Astronomy , McGraw-Hill Inc. 

Papoulis, A. (1965), Probability , Random Variables , and Stochastic Processes , Chapter 10, McGraw-Hill, Inc. 
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APPENDIX F 


THE SYNTHESIZED BEAM 

The synthesized beam G(q mn ) as defined in the text (Equation 24a) is: 

G(q miri - q) = DFT [ W exp(-j27r 3 -q) ] (FI) 

where DFT denotes a discrete Fourier transform in two dimensions on the variables and w is an arbitrary 

“weight” added here to account for biases in the measurements (e.g., differences in the effective time-bandwidth product 
or redundancies) or to permit beam shaping to control side lobes. Since the exponential is separable in 17 * the synthesized 
beam G(q mn - q) also separates and can be factored into the product of two functions: G(q mn - q) = G m (q m - q x )« 

G n (q n - q y ). Each of the G m n (q m n - q x y ) is a one-dimensional, discrete Fourier transform of the exponential 
exp(-j27M7 m n q x ). Using the following definition of the discrete Fourier transform 

DFT[ F ] = (1/M) ^ F(s) exp(j2ir [k/M]s) (F2) 

s = {M} 

each of the one-dimensional transforms can be written as follows: 

G m (q m - q x ) = (1/M) W s exp(j2ir[(m/M)s - v x qj) (F3a) 

s = {M} 

G n (q„ - q y ) = ( l/N) ^ Y s exp(j2ir[(n/N)s - ij y q y ]) (F3b) 

s={N} 

Equations F3 can be evaluated explicitly in the important, special case when the measurements are made at (or can be 
mapped onto) baselines which are multiples of some minimum spacing d x and d y in each dimension. The sums in Equa- 
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tions F3 can be evaluated in closed form in this case and the results will be given below. However, since the sums are 
identical, except for an exchange of subscripts, the discussion below will be for only one term, G m (q m - q x ). 

When the measurements are made at multiples of some minimum spacing, d x , the parameter i? x (s) becomes an in- 
teger multiple of (d x /c)v 0 = d x /X (e.g., see Equation Bll, Appendix B). In this case, Equations F3 have the form: 


G m (q m - q x ) = (1/M) W s exp(jk 0 dJmX/(Md x ) - q x ]s) (F4) 

s = {M} 

where k 0 = 27 r/X. Comparing the right- and left-hand sides, it is clear that q m - mX/(Md x ) where M is the number of 
terms in the sum. 

In order to proceed further, consider the special case when all of the weights, W s , are equal. In this case, the sum 
in Equation F4 can be evaluated explicitly (it is a geometrical progression) and one obtains: 

1 sin[7tMd x (q m - q x )/X] 

Gm^m ~ <U> = “ ( F5 > 

M sin[7rd x (q m - q x )/X] 

sinc[Md x (q m - q x )/X] 
sinc[d x (q m - q x )/X] 

where sinc(x) = sin( 7 rx)/ x. Notice that Equation F5 is just the array factor of a linear array of M identical antenna ele- 
ments with uniform spacing, d x (Jordan and Balmain, 1968). 

Notice that there are M terms in the sum in Equations F3 and F4 and, therefore, also M possible values of q m . The 
M terms consist of measurements at M baselines, the complex conjugate of these measurements and a measurement at the 
zero baseline. Hence, one has M = 2 M + 1 terms in the summation in Equation F4. The values of the index m are: m 

A A 

- 0, ±1, ±2 ... ±M. Also, notice that the physical length of the synthesis array (i.e., the maximum baseline) is Md x ; 
but the equivalent electrical length as compared to a linear antenna array with the same resolution is Md x . 

Also, notice that each value of m will result in a data point (pixel) in the image (M points in the x-dimension and N 
points in the y-dimension). The increment between these data points in direction cosine space is q m . However, these 
points are not separated by a full synthesized beam width. To see this, define the resolution of the synthesized beam to be 
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the distance between the first nulls on each side of the main lobe. Denoting this distance by q,,, it is clear from Equation 
F5 that: 

q a = 2\/ Md x = X / D x ( F6 ) 

where D x = M d x is the physical length of the array used to synthesize the beam (i.e., the maximum baseline). Since -1 
< q x < + 1 it follows that there are 2/q^ = Md x /X = 2D x /X contiguous resolution cells in the image. On the other hand, 
there are a total of M pixels in the reconstructed image in each dimension and each pixel is separated by q m = X/(Md x ) 
(Equation F4). Since q m = q c /2 it follows that there are two pixels per resolution cell. 

Finally, notice that with uniform weight, grating lobes will appear in the image at integer multiples of X / d x = q x 
(Equation F5). Since — 1 < q x < + 1, one must have X > d x in order to avoid the presence of grating lobes in the 
field of view of the sensor. 

REFERENCES: 

Jordan, E.C. and K.G. Balmain (1968), Electromagnetic Waves and Radiating Systems, Chapter 12, Prentice-Hall, 
Englewood Cliffs, New Jersey. 
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